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(CDICI) as described by Drummond and Rosenbluth (Ref. 

6 ). 


ABSTRACT 

Broadband electrostatic turbulence has been observed re- 
cently in the auroral regions in conjunction with transverse 
d.c. electric fields. Generally, field-aligned electron drifts are 
also found in this region. Using a spatially two-dimensional 
electrostatic PIC simulation code, we have studied the sta- 
bility of a plasma equilibrium characterized by a localized 
transverse d.c. electric field and and a field-aligned drift for 
L <C T x , where L x is the simulation length in the x direction 
and L is the scale length associated with the d.c. electric 
field. It is found that the d.c. electric field and the field- 
aligned current can together play a synergistic role to enable 
the excitation of electrostatic waves even when the thresh- 
old values of the field aligned drift and the E x B drift are 
individually sub critical. The simulation results show that 
the growing ion waves are associated with small vortices in 
the linear stage, which evolve to the nonlinear stage domi- 
nated by larger vortices with lower frequencies. 

Keywords: Parallel Currents, Transverse Electric Fields, In- 
homogeneous Energy Density Driven Instability, Vortices, 
E x B drift, Anomalous Diffusion. 


1. INTRODUCTION 

Recently, shears in the transverse flow velocity and associ- 
ated broadband electrostatic turbulence have been observed 
in a wide range of the auroral region (Ref. 1-5). Basu et ah 
(Ref. 2,3) using DE 2 satellite data and Earle (Ref. 4) using 
COPE rocket data provide clear evidence of strong velocity 
shears (10 to 25Hz) and associated broadband fluctuations 
in the auroral F region. Besides the shear in the transverse 
flow velocity, a magnetic field-aligned current (FAC) also 

exists in this region. Generally the observed magnitudes 
of the FAC are below the threshold value necessary for the 
excitation of the Current Driven Ion Cyclotron Instability 
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Ganguli et al. have developed a kinetic theory (Ref. 7,8) 
to explain the electrostatic waves associated with transverse 
velocity shears. Recently the initial kinetic theory (Ref. 7,8) 
has been generalized to include continuous gradients in the 
d.c. electric field. In the weak shear limit (i.e V E (x)/Hi < 1, 
where V' E is the spatial derivative of the E X B drift and Hi 
is the ion cyclotron frequency), the distinctions between the 
Inhomogeneous Energy Density Driven Instability (IEDDI) 
and the kinetic K-H instability have been studied. The 
growth rates of the kinetic K-H modes were obtained for 
various values of u = k\\/k v where A;|| and k y are the compo- 
nents of the wave vector in the z (parallel to the magnetic 
field B) and the y directions. It was found that the K-H 
mode is localized below b = k*p] — 1 and k y L ~ 1, where L 
is the scale length associated with the d.c. electric field and 
Pi = Vi /Hi is the ion gyro- radius. Furthermore, the growth 
rates are reduced strongly as u increases. The growth rates 
for the IEDDI were also obtained. This mode becomes un- 
stable where b > 1 and k y L > 1. The unstable domains for 
both modes are distinct. 

Nishikawa et al. (Ref. 10) investigated the IEDDI by means 
of particle simulations. The linear analysis of waves ob- 
tained from the simulations show that the mode (0, 4) which 
corresponds to b = 3.38 and k y L = 4.31 has the maximum 
growth rate and falls on the IEDDI branch. The emphasis 
in this simulation (Ref. 10) was to demonstrate the exis- 
tence of the new mode (IEDDI) and distinguish it from the 
well-known Kelvin- Helmholtz branch and the kinetic modi- 
fication of this branch. 

Subsequently, Ganguli and Palmadesso (Ref. 11) and Gan- 
guli et al. (REF. 12) have shown that the simultaneous exis- 
tence of transverse localized electric (TLE) fields and field- 
aligned currents (FAC) can excite electrostatic waves even 
though the electric field and the current are individually sub- 
critical. In the presence of TLE fields, the Doppler shift due 
to the E x B drift plays an important role and the thresh- 
old value of the FAC for exciting the CDICI is effectively 
lowered (Ref. 12). The character of the waves depend on 
R (= kyVE/k t Vd) y the ratio of the Doppler shifts due to the 
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drifts originating from the TLE and the FAC. For R < 1 
the waves are CDICI-like while for R > 1 the waves are 
IEDDI-like. 

Here we describe the numerical simulations that have been 
performed in order to understand the excitation of electro* 
static waves in the simultaneous presence of a field-aligned 
electron drift and a transverse d.c. electric field (Ref. 12) 
and it’s evolution to the nonlinear regime. In Sec. II. the 
simulation model and results are presented and concluding 
remarks and discussions are described in Sec. III. 


2. SIMULATION MODEL AND RESULTS 


A two-dimensional electrostatic code is used which retains 
the full dynamics of the ions in three dimensional velocity 
space. Electrons are treated by the guiding center approxi- 
mation in the perpendicular direction while the parallel mo- 
tion is treated exactly. We use a system length L x = 128A, 
L y = 32A, where A is the grid size which is equal to the 
electron Debye length, A c and n c A 2 = 36. In the previous 
work (Ref. 10), we found that the IEDDI has the max- 
imum growth rate around k y pi = 1.84 which corresponds 
to (0, 2) mode in the case of L y = 32A. By using the 
system length L y = 32 A instead of 64 A, we can isolate 
only the IEDDI (Ref. 10) by cutting out the KH instability 
from the system. The external electric field is applied in 
the form of E ox (x) = £’ <>r sech 2 [(x - 64)/L v ] in the x direc- 
tion which produces E x B drift in the y direction given by 
Ve{x) ~ — cE OT (x)/ B a . The parallel electron drifts are ini- 
tialized in the form of Vd(x) = Vgsech 2 [(;r— 64)/£ c ]. Here we 
maintain L c = L v , but this restriction can easily be relaxed. 


For the present problem, in which a nonuniform d.c. elec- 
tric field transverse to the uniform magnetic field is present, 
an appropriate initial distribution function has been con- 
structed out of the constants of motion (Ref. 9). For the 
purpose of initial loading in a computer simulation, the dis- 
tribution can be expressed in terms of the physical position 
x (for details see Ganguli et al. (Ref. 9)), 


2?r/o ± = 




[1 + 0(e)]. (1) 


Here tj = l^Vg(x)/Qi and v t = is the thermal velocity. 

Note that the distortion in the gyro-orbit introduced by the 
sheared d.c. electric field leads to a sustainable difference 
in the temperature in the two dimensions transverse to the 
uniform magnetic field. For weak shears (i.e. | V^/B |< 1), 
Eq. (1) reduces to, 


2*/oi = /?{1 + ^p[(f v “ v e(*)) 2 " ”«]} 

* ex P {“f [*£ + ( v v ~ ^(z)) 2 ]}- 


We see that for weak shears the use of a shifted Maxwellian 
for the initial loading is acceptable and the system will ad- 
just to the above distribution automatically with only minor 
relaxation (Ref. 10). If the Vg/Q correction is ignored, then 
Eq. (2) becomes identical to a Maxwellian shifted by the 


magnitude of the E x B drift in the y component of the 
velocity. This describes the distribution function transverse 
to the magnetic field. We use / 0 || = (^/2x) 1 / 2 exp (— — 
Vj(x)) 2 /2) as the distribution function along the magnetic 
field so that /<, = /oi/oj|- 

We first choose three set of parameters: (i) Vf = — 0.4u c 
and Vg = 0.0, (ii) Vf = 0.0 and Vg = 0.15?;,, and (iii) 
Vf = — 0.4t; e and Vg = 0.15v,-, where VJ, Vg, v e and Vi 
are the peak values of V^(x) and Vg(x), and the electron 
and ion thermal velocities, respectively. Other parameters 
used are B oy /B 0 (= k^/k y ) = 0.005, = 4, Ti/T e = 

3.5, m t /m e = 100, a = E% x /4irn e T e = 0.0125, e — pi/L = 
0.42, and L = L v = L c = 11 A. We display in Figure 1 
the time evolutions of the fluctuating electrostatic potentials 
and their corresponding power spectra for the mode (0, 2). 
The mode (0,2) corresponds to b= 3.38, k y pi = 1.84, k y L = 
4.31. Here the magnitude of the peak shear Vg/Ltti = 0.063. 
Also, since R = 14 >> 1, the waves are IEDDI-like. 



Figure 1. The comparison of the time evolutions and 
their spectra of the (0, 2) modes among the 
three cases. 

(a) and (b): Vf = -0.4v e and Vg = 0.0, 

(c) and (d): Vf = 0.0 and Vg = 0.15?;,-, 

(e) and (f): Vf = -0.4v e and Vg = 0.15^. 


As shown in Figure 1, the cases with only the parallel elec- 
tron drift (a) or only the transverse electric field (c), can- 
not sustain any growing wave. In their spectra (Figure 1 
(b) and (d)), the thermal noises at the cyclotron frequency 
and its higher harmonics are observed at low intensity. On 
the other hand, we see that the combination of the TLE 
and the FAC can sustain growing waves. The waves be- 
gin to grow and emerge from the background thermal noise 
around = 220 and subsequently reach the nonlinear 
stage around B,-/ = 300. As shown in Figure 1 (f), the 
peak of the frequency spectrum is located around —0.20,. 
In Figure 1 (f) we also find a number of smaller peaks around 
higher harmonics. 

In Figure 2 we show the time evolutions of the modes (0,3) 
and (0,4) ((a) and (c)) and their corresponding power spec- 
tra ((b) and (d)). Here rest of the parameters are unchanged 
from the previous case. 
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Figure 2. The time evolutions ((a), (c)) and their power 
spectra ((b), (d)) of perturbations of 
electrostatic potential for two modes with 
the parallel currents and the electric fields. 

(a) and (b): fc||A e || = 0.005 x 2rx 3/32 = 2.95 x 10" 3 , 
Jfc v A,|| =2rx 3/32 = 0.589. 

(c) and (d): Jfc||A e || = 0.005 x 2* x 4/32 = 3.93 x 10" 3 , 
fc„A e || = 2?r x 4/32 = 0.785. 


8- 887-82* 


<o) Ity*l60 



Figure 4. The time evolution of the contour plot of the 
of the electrostatic potential in the linear stage 
flit = 160 and 240, and in the nonlinear stage 
flit = 320, 328, and 336. 


The mode (0, 3) has a peak with the larger real frequency 
at -0.3H, as shown in Figure 2 (b). It is interesting to note 
that the power spectra (b and d) are remarkably similar 
in general features (e.g., spiky nature, generation of higher 
harmonics etc.) to the high time resolved spectra given by 
Prikryl et al. (Ref. 13), which were obtained from their 
radar measurement of the plasma irregularities in the auro- 
ral region. 


*-087-827-1 



Figure 3. The average ion flow velocity r y (x) (a) and 
the parallel electron currents V& (b) at 
flit = 0 and 320. 


The average ion flow velocities V y (x) and average electron 
flow velocities Vd(x) are shown in Figure 3. In the linear 
stage (f lit < 100) there is hardly any reduction in the peak 
of the velocity profile or any relaxation in the topology of 
the profile. At flit — 320 (thicker curves), the peaks are 
moved at x = 80 due to nonlinear effects. It should be 
noted that the value of peak of electron currents remains al- 
most unchanged. On the other hand, in the late quasilinear 
(flit ~ 240) and the nonlinear (flit > 320) stages the wave 
amplitude become large and can interact with and reduce 
the magnitude of the d.c. electric field. In the final stages 
of the simulation we see a reduction in the magnitude of 
Ve(x) : Since the real frequency is proportional to k y VE(x) 
(see Ganguli et al. (Ref. 9)), we find that in the quasilin- 
ear stage the real frequency starts reducing along with the 
magnitude of Vje(x). Further reduction in the real frequency 
occurs in the fully nonlinear phase. 

The evolution of vortices at the linear and nonlinear stages 
(flit > 160) is shown in Figure 4. The electrostatic potential 
in the x-y plane has strong vortices at the center at flit = 
240. These vortices move along the direction of the E X B 
drift. As shown in Figure 4, the excited waves are well 
localized around the center. In the later nonlinear stage the 
number of the vortices is increased and one of the vortices 
migrates outward from the center. 

The diffusion of the guiding center of ions in the case of 
the combination of TLE and FAC is shown at flit = 0, 80, 
160, 240, and 320 in Figure 5. A subset of ions with initial 
positions such that 58 < x < 64, are followed in time, and 
the local guiding centers of ions are plotted at each time. 
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Figure 5. The diffusion of the gyro-center of ions 

whose initial actual positions are 58 < x < 64 
at = 0, 80, 160, 240, and 320. 


At the time (ft = 160) when the waves begin to grow 
out of the thermal noise, slight waving of the column is 
observed. At the beginning of the nonlinear stage Qit ~ 240, 
the column is twisted. In the nonlinear stage we see strong 
diffusion of the ions. 

The local electric fields E ZtV at x = 74, y = 14 is diagnosed, 
and the power spectrum is obtained. 


A-G87-83I-I 



Figure 7. The average field energy ( E*) t and the d.c. 

type polarized electric field energy -(£*)? are 
localized around the center of the system. 


The profiles of overall wave energy are shown in Figure 7. 
The time-averaged excited wave energy (Ef) t is localized 
eu-ound the center of the system. Obviously the excited wave 
energy is much larger than the thermal energy outside the 
electric field region. The time- averaged polarization electric 
field energy with minus sign ~{E X )1 is also plotted, which 
is much smaller than the excited wave energy. 

3. DISCUSSIONS 
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Figure 6. The time evolution of local electric field E XyV 
at x = 74, y = 14 (a) and its power 
spectrum (b). 


We have investigated the electrostatic waves driven by the 
combined effects of a localized transverse electric field and 
parallel electron drifts by means of simulation with the as- 
sistance of the nonlocal kinetic theory (Ref. 12). 

We have performed a number of simulations for this insta- 
bility. Simulation results show that electrostatic waves are 
excited in the regions where the ExB drift is localized 
in the simultaneous presence of parallel electron drifts and 
transverse electric fields as predicted (Ref. 12). The mag- 
nitudes of the FAC and the TLE were chosen to be indi- 
vidually subcritical so that simulations with the FAC or the 
TLE when considered in isolation shows no instability, but 
when the combination is considered we find clear signatures 
of a growing wave. The physical mechanism underlying this 
instability is elucidated in the companion paper (Rerf. 12). 

The nonlinear phenomena such as diffusion and coalescence 
of vortices are investigated. In the linear stage, smaller vor- 
tices are generated and larger vortices with the lower fre- 
quencies are dominant in the nonlinear stage. In the non- 
linear stage the ions diffuse strongly due to the large scale 
vortices. 


As shown in Figure 6 (a), around Hit = 230 sharp peaks are 
produced. The spectrum has several peaks around 0.40,- 
which correspond to the real frequencies in the linear stage. 
The largest peak is produced by the vortices in the nonlin- 
ear stage. 


Our simulation results may explain the excitation of electro- 
static waves of the observations (Ref. 2-4). However, more 

simulations are necessary to explain various other features 
of the observations. We are currently improving our simula- 
tion code to include the neutral-species collisions which are 
essential for the lower auroral region studies. 
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